Proton-proton correlations observed in two-proton decay of 19 Mg and 16 Ne 
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Proton-proton correlations were observed for the two-proton decays of the ground states of 19 Mg 
and 16 Ne. The trajectories of the respective decay products, 17 Ne+p+p and 14 0+p+p, were mea- 
sured by using a tracking technique with microstrip detectors. These data were used to reconstruct 
the angular correlations of fragments projected on planes transverse to the precursor momenta. 
The measured three-particle correlations reflect a genuine three-body decay mechanism and allowed 
us to obtain spectroscopic information on the precursors with valence protons in the sd shell. 

PACS numbers: 21.10.-k; 21.45.+v; 23.50.+Z 
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The recently discovered two-proton (2p) radioactivity 
is a specific type of genuine three-particle nuclear decays. 
It occurs when a resonance in any pair of fragments is lo- 
cated at higher energies than in the initial three-body 
(p+p+"core") nucleus, and thus simultaneous emission 
of two protons is the only decay channel. Three-body 
systems have more degrees of freedom in comparison with 
two-body systems, hence additional observables appear. 
In the case of 2p emission, the energy spectra of sin- 
gle protons become continuous, and proton-proton (p-p) 
correlations are available, which makes them a prospec- 
tive probe for nuclear structure or/and the decay mecha- 
nism. For example, the first p-p correlations observed in 
the 2p radioactivity of 94 ™ l Ag have revealed strong pro- 
ton yields either in the same or opposite directions which 
called for a theory of 2p emission from deformed nuclei 
[l| . Two-proton emission can also occur from short-lived 
nuclear resonances or excited states (see, e.g., @, H, 0])- 
Though in this case the mechanism of 2p emission may 
depend on the reaction populating the parent state, such 
nuclei can be easily studied in-flight. E.g., the cases 
of 6 Be d [B| and 16 Ne Qwere studied by analyzing their 
p-p correlations in the framework of a three-body partial- 
wave analysis developed for three-particle decays of light 
nuclei. In particular, the study of 6 Be revealed the ex- 
istence of three-particle p+p+cv correlations [2j which 
matched the three-body components found theoretically 
in the p-shell structure of 6 Be Q . Very recently, p-p cor- 
relations were also observed in 2p radioactivity of 45 Fe 
[1 Q where both the lifetime and p-p correlations were 
found to reflect the structure of p/-shell 2p precursors @ . 
Such a way of obtaining spectroscopic information is a 
novel feature compared to studies of two-particle decays. 



In the present paper, we study for the first time the 
p-p correlations in sd shell nuclei via examples of the 2p 
decays of 19 Mg and 16 Ne. These nuclei with very different 
half lives (T 1/2 « 4-KT 9 s [To| and T 1/2 « 4-KT 19 s 
[llT ] , respectively) and presumably different spectroscopic 
properties may serve as reference cases illuminating the 
nuclear structure of other possible 2p emitters with sd- 
wave configuration. 

The decay properties of the 16 Ne and 19 Mg ground 
states and the related resonances in 15 F and 18 Na are 
shown in Fig. [1] which compiles the data from Refs. [H 
EH, El El EL IS El| and this work. The ground states 
of both isotopes decay only by simultaneous 2p emis- 
sion while their excited states are open for sequential lp 
decays via intermediate unbound states in 15 F and 18 Na. 

The quantum-mechanical theory of 2p radioactivity 
which uses a three-body model [IT], EH Ell, predicts 
the p-p correlations to be strongly influenced by nu- 
clear structure together with Coulomb and three-body 
centrifugal barriers. In particular, the newly discovered 
2p-radioactivity of 19 Mg [1(1 was predicted to be charac- 
terized by p-p correlations which reflect the sd configura- 
tions of the valence protons [2(1 . A similar effect is found 
in 1 Ne, where the s-wave configuration was predicted 
to dominate, contrary to its mirror 16 C, thus breaking 
isospin symmetry [2l|. A complementary approach 
in describing 2p decays is the mechanism of sequential 
emission of protons via an intermediate state (see e.g., 
[22l]). It includes also the traditional quasi-classical di- 
proton model with emission of a 2 He cluster, assuming 
extremely strong p-p correlations [23l |24| . The predic- 
tions of these models differ dramatically with respect 
to the p-p correlations, suggesting them as a sensitive 
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FIG. 1: States observed in 16 Ne, 19 Mg and the corresponding 
intermediate systems 15 F, 18 Na. Decay energies (in keV) are 
given relative to the respective p and 2p thresholds. Most 
values have been taken from the literature (Refs. [ToL [Til . Il2l . 
[l3.ll5l.ll6l]). those in bold print are from the present work. 
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FIG. 2: Left panel: Cartoon of momentum correlations 
kpi-Hi — k P 2-Hi expected for a direct three-body decay (the 
grey area labeled A;3_r) and sequential 2p decay (grey boxes 
with black peaks labeled k2n). Upper-right panel: A sketch 
of the kinematical enhancement of an angular p-HI correla- 
tion at the maximum possible angle for a given momentum 
between the decay products. Lower-right panel: the corre- 
sponding angular p-HI distribution. 



probe of the 2p-decay mechanism (see the detailed pre- 
dictions below). 

Our experiment to investigate 2p-emission from 19 Mg 
and 16 Ne was performed by using a 591^4 MeV beam of 
24 Mg accelerated by the SIS facility at GSI, Darmstadt. 
The radioactive beams of 20 Mg and 17 Ne were produced 
at the Projectile-Fragment Separator (FRS) [25[ with av- 
erage intensities of 400 and 800 ions s" 1 and energies of 
450^4 and 410^4 MeV, respectively. The secondary 1- 
n-removal reactions ( 20 Mg, 19 Mg) and ( 17 Ne, 16 Ne) oc- 
curred at the mid-plane of FRS in a secondary 2 g/cm 2 
9 Be target. Special magnetic optics settings were ap- 
plied, the first FRS half being tuned in an achromatic 
mode using a wedge-shaped degrader, while its second 
half was set for identification of the heavy ions (HI) with 
high acceptance in angle and momentum. 

A sketch of the experimental set-up at the FRS mid- 
plane has been shown in Fig.l of Ref. [l(| and explained 
in detail there. A microstrip detector array [26| consist- 
ing of 4 large-area (7x4 cm 2 ) double-sided silicon detec- 
tors (with a pitch of 0.1 mm) was positioned downstream 
of the secondary target. This array was used to measure 
energy loss and position of coincident hits of two protons 
and a heavy fragment, thus allowing us to reconstruct 
all decay-product trajectories and derive the coordinates 
of the reaction vertex and the angular p-p and proton-HI 
correlations. The conditions to select the true HI+p+p 
events were: (i) a minimal distance between the proton 
and heavy ion trajectories of less than 150 /im, and (ii) a 
difference between the two longitudinal coordinates of the 
vertices defined by two p-HI pairs (taken from the same 
HI+p+p event) within the range defined by the experi- 
mental uncertainty of 0.3-1 mm depending on detection 
angle. The achieved angular resolution in tracking the 
fragments was ~1 mrad. More details concerning the de- 
tect or p erformance and tracking procedure are given in 
floL l27l . [2R [29i | . Another position-sensitive silicon de- 
tector and a multi-wire chamber were used upstream of 



the target for tracking the 20 Mg( 17 Ne) projectiles. The 
heavy 2p-out residuals ( 17 Ne and 14 0) were unambigu- 
ously identified by their time of flight, magnetic rigid- 
ity, and energy loss measured with the position-sensitive 
scintillator detectors at the second half of FRS. 

The 2p decays of the ground states of 19 Mg and 16 Ne 
were identified by using angular correlations between the 
single protons and their respective cores 17 Ne and 14 
which allowed measuring the 2p-decay energies. This is 
analogous to identifying a reaction via the scatter plot 
illustrated in Fig. [5] (left panel). In general, 2p-decay 
may proceed via either sequential or direct decay mech- 
anisms. The first case can be described as two consecu- 
tive lp decays, with the p-HI spectra reflecting the re- 
spective p-HI resonances [12]. Their (pi-HI)-(p2-HI) 
scatter plot should display the sequential and direct de- 
cay events in the respective kinematical areas marked in 
Fig. d]by the relative momenta k2R- In the second case, 
the simultaneously emitted protons are likely to share 
the 2p-decay energy evenly, with both p-HI spectra be- 
ing identical and peaked at E/2 [l7], 0|- In this case, 
the area marked in Fig.f5]should be populated, along 
the arc area with the root-mean-squared proton momen- 
tum being constant. Sequential proton emission from a 
single 2p-parent state via narrow p-HI resonances should 
yield double peaks while 2p de-excitation of continuum 
parent states with p-HI final state interactions should 
reveal "slices" as shown in Fig. [5J 

Similar structures can be found in the angular p i-hi — 
6 P 2-Hi correlations due to the following reason. Because 
of a strong kinematic focusing at intermediate energies, 
the lp decay leads to a characteristic angular p-HI cor- 
relation when the proton ejected isotropically in the lp- 
precursor frame is emitted within a narrow cone around 
the HI with the maximum intensity around the largest 
possible angle (Fig. [2j right panel). The p-HI angles 
reflect the transverse proton momentum relative to the 
HI one, and are therefore correlated with the precur- 
sor's decay energy. Thus sequential 2p decays from ex- 
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FIG. 3: (a) Angular (pi- 0)-(p2- O) correlations obtained from the measured O+p+p events ( colored boxes with scale 
shown on right). The grey areas indicate simultaneous and sequential 2p emission in analogy to those shown in Fig. [2] (b) 
Angular p- 14 distribution (full circles with statistical uncertainties) obtained from the data shown in panel (a) by selecting 
the other proton angle 6 P 2-o within the range from to 45 mrad, which corresponds to the ground state of 16 Ne. The solid 
curve represents the Monte Carlo simulation of the detector response for 16 Ne 9 .s— > 14 0+p+p with a 2p-decay energy of 1.4(1) 
MeV [1^1 . The dashed line is a sum fit to the data. The dotted curve indicates the background as explained in the text, (c) 
Angular p- 14 distribution obtained from the data shown in (a) by selecting the other proton angle from 120 to 150 mrad, 
which corresponds to p- 14 final-state interactions due to the ground and excited states of 15 F. The dashed and dotted curves 
are the Monte Carlo simulations of the known lp decays of the ground and first excited states of F with the lp-decay energies 
of 1.56(13) and 2.85(4) MeV, respectively [l3j]. The dash-dotted and dash-dot-dotted curves show the result of fits to the 
p- 14 correlation, indicating unknown excited states in 15 F with lp-decay energies of 4.9(2) and 6.4(2) MeV, respectively. 
The solid line is the sum fit. 



cited states in parent nuclei should be mostly located 
in peaks with tails along the respective slices in the an- 
gular 6pi-Hi — 9p2-Hi correlations, in analogy to those 
sketched in Fig. [2] (left panel). In direct 2p decays, the 
single-proton energy spectrum always exhibits a rela- 
tively narrow peak centered close to half of the 2p-decay 
energy; such energy distribution is a stable feature of 
this decay mechanism [20l . [23| . Correspondingly, a bump 
should appear in the angular correlations in the same 
way as it should appear along the arc marked by k^R in 
the scatter plot in Fig. [2] This correspondence between 
angular and momentum correlations has been used to de- 
rive the 2p-decay energy of 19 Mg (see Fig. 2,4 in [l(| and 
the respective discussions). 

For the 2p-decay of Ne, the angular correlations of 
each coincident proton with respect to the 14 O momen- 
tum, p i-o — 6p2-o, derived from the measured 14 0+p+p 
coincidence events, are shown in Fig. [3]Ja). The events 
with the smallest angles fall into a distinct cluster around 
#p_o=35 mrad while most of the other events are located 
in the slices centered around 70 and 95 mrad. These two 
groups can be attributed to the direct 2p decay from 
the 16 Ne ground state and to the sequential emission of 
protons from excited states in 16 Ne via the 15 F ground- 
state, respectively. We shall refer to these events as the 
"ground state" and "excited state" events, respectively. 
The latter group includes also events resulting from the 
fragmentation reaction 17 Ne^ 14 0+p+p+n. 

To disentangle the "ground state" from the "excited 
state" events, we made a slice projection from the mea- 
sured (pi- 14 0)-(p2- 14 0) correlations in Fig.[3Ja) by se- 
lecting the angle of one of the protons within the range 
0-45 mrad, where the 2p decay of the 16 Ne ground state 



is expected to show up. Figure \5j[b) displays the angular 
correlations # p i_o corresponding to the "ground state" 
gate in the other pair # P 2-o- The peak around 35 mrad 
(the suggested "ground state" ) dominates the spectrum, 
whereas few correlations can be seen between a proton 
from the "ground state" and another proton at larger an- 
gles. This means that the two protons from the "ground- 
state" are correlated, i.e. this peak can be explained by 
an emission of protons from the ground state in 16 Ne. 

For a more quantitative analysis, the data are com- 
pared to a Monte Carlo simulation of the response of 
our setup to the direct 2p-decay 16 Ne— > 14 0+p+p with 
the known 2p-decay energy of 1.4(1) MeV [HJ by using 
the GEANT program [3(J. The simulations took into 
account the above-mentioned experimental accuracies in 
tracking the fragments and in reconstructing the ver- 
tices, trajectory angles etc. The normalized simulation 
reproduces the data in the low-angle peak very well. The 
contribution from the "tail" of the higher states to the 
ground-state peak amounts to about 20%. The shape 
of this distribution is assumed to have the same shape 
as the 9 p i-o distribution selected within the 9 P 2~o range 
just outside the ground-state region, from 48 to 160 mrad 
(the dotted curve). 

Figure \%lc) displays an example of an angular p- 14 
distribution for "excited states" obtained from Fig. [3ja) 
by selecting the angular range of the other proton from 
120 to 150 mrad, which corresponds to p- 14 final state 
interactions due to the ground and excited states of 15 F. 
The Monte Carlo simulation of known one-proton decays 
i5p_^i4Q_|_p £ g roun( j an( j fi rs t excited states in 15 F 

with the lp-decay energies of 1.56(13) and 2.85(4) MeV 
[HI], respectively, reproduces well the two lowest-angle 



4 




E p . p IE cos(6 ( .) or cos(9 1 ') 



FIG. 4: (a) Three-body correlations predicted by the three- 
body model 20] for the 2p-decay of 19 Mg, plotted as a func- 
tion of the relative energy between two protons, E p - P /E, and 
the angle 6^ between the relative momenta of the fragments as 
illustrated in the figure. Extreme cases of p-p correlations are 
sketched and related to the corresponding areas in the corre- 
lation plot, (b) The p-p energy spectra from the 2p-decay of 
19 Mg calculated for different weights W of s-p-d shell config- 
urations in 19 Mg. (c) Typical intensity distributions plotted 
as a function of 8k (see Fig. |4ja)) in the rest frame of the 
2p-precursor (solid curve) and its analog in the lab system 0' k 
projected on the transverse detector plane (dashed curve). 



peaks. The two higher-lying peaks indicate lp decays 
of unknown excited states in 15 F with derived lp-decay 
energies of 4.9(2) and 6.4(2) MeV. The excited states in 
19 Mg, 16 Ne, 18 Na and 15 F will be addressed elsewhere. 

We turn now to the discussion of angular p-p correla- 
tions following 2p decays. When the spin degrees of free- 
dom are neglected and the total decay energy E can be 
considered as fixed, the three-body correlations are com- 
pletely described by two variables. A convenient choice 
is an energy-distribution parameter E p ^ p /E (E p ^ p is 
relative energy between two protons) and an angle 9k 
between the relative momenta of the fragments. FigfJ] 
shows such distributions predicted by the three-body 
model for the 2p-decay of 19 Mg [H,^- The three-body 
model predicts a distinctive correlation pattern which 
features an enhancement at small E p _ p due to final-state 
interaction and a suppression in the regions of strong 
Coulomb repulsion (E p _ p /E ~ 0.5, cos(9k) ~ ±1). The 
predicted energy distributions are sensitive to the struc- 
ture of the precursor [Fig. HJb)]. Similar predictions are 
available for 16 Ne |2lj . 

In our experiment, we were able to measure the open- 
ing angle Q p - P between protons whose distribution re- 
flects the E p _ p correlations. FigfS] shows the exper- 
imental angular p-p distributions obtained from triple 




FIG. 5: Angular p-p correlations from the 2p-decays of the 
ground states of 19 Mg (left) and 16 Ne (right)) obtained from 
the measured 14 0+p+p and 17 Ne+p+p events, respectively, 
by selecting both protons from the respective p- 14 0( 17 Ne) 
angular ranges. The solid curves show the fits to the data, 
obtained by using the three-body model calculations (54% 
and 88% of d-wave configuration in 16 Ne and 19 Mg, respec- 
tively.) The dotted curves show the background contributions 
estimated as described in the text. The dashed curves are the 
diproton model predictions, and the dash-dotted curves are 
the phase-space simulations of the 2p-decays (isotropic pro- 
ton emission in the 2p-precursor's rest frame). 

14 0-r-p+p and 17 Ne-r-p+p events gated by the conditions 
that both protons originate from the "ground states" 
of 16 Ne and 19 Mg. These gates were inferred from the 
respective angular p-HI correlations as discussed above. 
The events representing the "ground state" 2p decay ac- 
tually contain contributions from the ground state and 
background contributions from both, excited states of the 
parent nucleus and fragmentation reactions. Therefore, 
using the angular-correlation data discussed above, we 
empirically evaluated the shapes of the background com- 
ponents by projecting triple events with the p-HI gates 
shifted away from the "ground state" region towards 
larger angles. The resulting p-p background contribu- 
tions shown by the dotted curves in Fig. [5] constitute 
about 20 % of all p-p correlation data for 16 Ne and 25 
% for 19 Mg (see Fig. EJb) and Fig. 4 (c) in [HI); they 
were subtracted from the original p-p correlations. As 
one can see in Fig. the predictions following from the 
assumption of a diproton emission fail to describe both 
the 16 Ne and 19 Mg data while the three-body model re- 
produces the shapes of both distributions. In the 19 Mg 
case, the best description is obtained with the <i-wave 
configuration dominating. The 16 Ne data give evidence 
for nearly equal s- and d-wave components. 

In Fig. [6l the intensity distributions are displayed as 
a function of cos(^). The angle 9' k (see Fig. H£a)) was 
defined by a line connecting the two points where two 
protons hit the same detector and by a vector join- 
ing the middle between the 2p hits and the point of 
a related heavy-ion hit, in analogy with the angle 9 k 
shown in Fig. [3Ja). The typical theoretical prediction 
for such a distribution is shown in Fig. [3Jc). The dipro- 
ton model predicts flat angular distributions in contrast 
to the experimental data in both cases. Only the three- 
body model can reproduce the characteristic shapes of 
the observed correlations with the broad bumps around 
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FIG. 6: Three-body correlations from the 2p-decays of the 
ground states of 16 Ne and 19 Mg (full circles with statisti- 
cal uncertainties) obtained from the measured 14 0+p-|-p and 
17 Ne+p-|-p events by selecting an angle 9' k (see text). The 
solid curves are the three-body model calculations (assuming 
54% and 88% of d-wave configuration in 16 Ne and 19 Mg, re- 
spectively) . The dashed curves are the diproton model predic- 
tions, and the dash-dotted curves are the phase-space simula- 
tions illustrating the response to an isotropic 2p-emission in 
the precursor's rest frames (both simulations are not normal- 
ized). The dotted curves show the background contributions 
estimated from all measured HI+p+p events. 

cos{9' k )—Q (the indicated spikes at cos(^)~±l predicted 
by all calculations are less conclusive). Such a shape is 
a manifestation of the "Coulomb focusing" efficiently re- 
pulsing the fragments from large regions of the momen- 
tum space (see Fig.[4ta)). These distributions are weakly 
sensitive to the assumed structure of the parent states 
but are an exclusive feature of the three-body model. 

In summary, the measured three-particle correlations 
from the 2p decay of the ground states of 16 Ne and 19 Mg 



are described quantitatively by the predictions of the 
tree-body model [20J, in contrast to the quasi-classical 
"diproton" model which fails to describe our observa- 
tions. These correlations are sensitive to the structure of 
the decaying nucleus. Thus the comparison between ex- 
periment and theory allows one to obtain spectroscopic 
information about the parent states. In 16 Ne, the data 
are consistent with strong s/d mixing [2l[. In 19 Mg, 
the dominating ci-shell configuration is the preferable de- 
scription which is also consistent with the lifetime in- 
formation 10]. The method of measuring 2p-decays in 
flight by precisely tracking all fragments with microstrip 
detectors provides new specific observables, thus yield- 
ing valuable spectroscopic information on such exotic 
isotopes. Information about two-body subsystems, e.g., 
15 F, is simultaneously obtained. Systematic studies of 
other 2p emitters predicted theoretically [lj| [3l[ are fore- 
seen with this novel technique. 
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